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;  Abstract 

rirst  we  present  approximate  analytical  solutions  to  the  wave  equation  inside  an  over- 
moded  metallic  sectangular  waveguide.  The  cold  eigenmodes  are  expressed  in  terms  of 
cylindrical  Gauss ian-Hermite  functions  times  trigonometric  functions  to  insure  the  boundary 
conditions.  Next,  we  discuss  a  numerical  three-dimensional  calculation  for  a  Free  Electron 
laser  (FED  amplifier  which  is  based  on  the  Lienard-Hiechert  solution  of  the  .Maxwell's 
equations  cast  in  an  integral  form.  This  approach  is  readily  and  efficiently  extended  to 
include  the  effects  of  the  metallic  boundaries  of  the  waveguide  by  means  of  the  method  of 
"image  currents".  Finally,  the  radiation  field  in  the  cavity  emitted  by  the  electrons  in 
the  presence  of  t.he  combined  fields  of  a  co-propagating  eigenmode  wave  plus  a  linearly 
polarized  magnetic  undulator  is  expanded  in  terms  of  cavity  eigenmodes.  This  expansion 
allows  us  to  compute  the  gain  per  resonator  mode. 

Introduction  ' 


A  rather  extensive  literature  on  the  3-d  FEL  theory  is  presently  available.  Progress 
in  t.his  area  has  been  occurring  in  the  last  two  years  as  a  result  of  efforts  of  several 
groups  which  have  carried  out  numerical  as  well  as  semi-analytical  calculations ‘ ^ “ 

The  mam  objectives  of  these  calculations  are:  1)  to  incorporate  into  the  theory  the 
finite  transverse  dimensions  of  the  optical  field  so  to  be  able  to  handle  the  outstanding 
problem  of  the  filling  factor  that  plagued  the  more  naive  1-d  theories  2)  to  describe  the 
optical  quality,  i.e.,  angular  distribution,  frequency  content  3)  and  lastly  to  provide 
secure  ground  for  the  optical-resonator  design  for  FEL. 

In  this  paper  we  focus  our  attention  on  the  waveguide  resonator  of  the  FEL  experiment  at 
the  University  of  California  Santa  Barbara  (UCSB) .  Its  mechanical  design  and  main  character¬ 
istics  are  discussed  in  the  accompanying  paper. ^ 

Practical  resonators  for  FEL  must  contain  the  following  features  in  addition  of  providing 
t.he  containment  needed  to  propagate  the  electron  beam  along  the  undulator:  1)  a  means  to 
inject  and  extract  the  high  current  relativistic  electron  beam  2)  possibility  of  allowing 
relatively  large  amplitude  for  the  periodic  magnetic  field  (reducing  the  gap)  to  insure 
above- losses  small  signal  gain  3)  small  beam  waist  to  increase  the  small-signal  gam  of  the 
laser  4)  low- losses  resonator. 

Some  of  these  requirements  are  not  easily  met  by  conventional  open  resonators  in  the  sub- 
millimeter  region  of  the  spectrum.  An  appealing  alternative  is  to  use  waveguide  resonators° 
which  offer  the  advantage  of  small  undulator  gap  and  small  optical  mode  area. 

In  section  II. 1  we  discuss  a  metallic  waveguide  resonator'  which  meets  the  design  criteria 
numbered  above.  We  present  analytical  approximate  solutions  of  the  wave  equation  with 
metallic  boundary  conditions  corresponding  to  low-order,  overmoded,  rectangular  waveguide 
with  infinite  aperture  cylindrical  mirrors.  The  solutions,  eigenmodes  of  the  "cold* 
resonator  are  expressed  naturally  in  terms  of  cylindrical  Gaussian-Hermite  functions.  In 
t.he  following  section  II. 2  proceeding  according  to  the  standard  technique  of  first  assuming 
perfectly  conducting  walls,  we  calculate  the  waveguide  losses  from  the  tangential  component 
of  t.he  magnetic  field  at  the  walls.  We  show  that  the  n-modes  propagate  along  the  z-axis 
with  negligible  straight- losses  in  the  far-infrared  region. 

The  section  111  is  dedicated  to  t.he  numerical  calculation  of  t.he  radiation  field  produced 
oy  a  relativistic  electron  beam  inside  a  waveguide,  moving  in  the  combined  fields  of  a 
linearly  polarized  magnetic  undulator  and  a  propagating  eigenmode  of  the  cavity. 

All  3-d  FEL  theories  published ^ “  so  far  assumed  open  resonators,  i.e.,  the  presence 
of  the  metallic  walls  along  t.he  interaction  region  have  been  ignored.  However,  to  achieve 
.arger  undulator  magnetic  fields  it  is  imperative  to  reduce  t.he  undulator  gap  and  conse- 
guently  the  .-.eight  of  t.he  cavity. 

The  radiation  from  a  relativistic  electron  beam  wiggling  in  t.he  interaction  region  is 
concentrated  along  a  narrow  cone  in  t.he  axiai-direction  with  an  angular  aperture*  1.  '•'’TT. 


When  the  transverse  diinension  of  the  radiation  cone  is  comparable  to  the  height  of  the 
resonator  channel  then  additional  interferences  occur  due  to  the  reflections  at  the  walls, 
increasing  the  power  density  at  the  observat:.on  plane. 


To  account  for  these  reflections,  i.e.,  to  satisfy  the  metallic  boundary  conditions,  we 
turn  to  the  "image  current*  approach^  (rectangular  waveguide)  which  allows  an  efficient 
and  fast  way  of  computing  the  field  inside  the  cavity  from  the  radiation  of  the  electron 
beam  in  free  space. 

The  3-d  electromagnetic  fields  proxiuce  by  an  electron  beam  in  an  undulator  can  be  consid¬ 
ered  as  the  coherent  superposition  of  radiation  pulses  produced  by  each  individual  electron. 
Each  pulse  contains  the  information  about  velocity  and  acceleration  of  the  corresponding 
electron  and  is  analytically  represented  by  the  Lienard-Wiechert  solutions. - 

After  adding  these  solutions  and  extracting  the  slowly  varying  amplitude  and  phase,  by 
means  of  the  standard  averaging  procedure  over  many  wavelengths  of  light  we  cast  the  total 
radiation  field  as  an  integral  over  the  "electron*  time  expended  inside  the  undulator.  To 
solve  it  numerically  it  is  necessary  to  Icnow  the  dynamics  of  the  electrons.  This  is  pro¬ 
vided  by  the  pendulum  equation. 

In  section  III. 3  we  show  how  to  expand  the  total  radiation  field  inside  the  cavity  in  the 
Fresnel  approximation  as  a  superposition  of  the  "cold”  eigenmodes  of  the  resonator. 

In  the  last  section  III. 4  we  define  the  power  gain  per  mode  for  the  UCSB  waveguide 
resonator  and  we  plot  the  gain  ciirves  as  function  of  the  resonance  parameter  for  the  first 
t.hree  Gaussian  transverse  modes.  We  discuss  the  relative  amplitude  of  the  maximum  of  the 
gain  curve  and  its  position. 


II.  Cold  eigenmodes  of  a  waveguide  resonator 

1.  Analytical  approximate  solutions  of  the  wave  equation  with  boundary  conditions 

The  UCSB  FEL  waveguide  is  a  metallic  rectangular  cross  section  channel  that  allows  both 
small  undulator  gap  and  small  radiation  field  beam  waist.  The  geometry  of  the  problem  is 
shown  in  Figure  I  as  well  as  the  approximate  dimensions  of  the  cavity. 


Figure  1.  A  schematic  drawing  of  the  UCSB 
FEL  resonator.  a»13.34  cm, 
b«1.9  cm,  d«10.0  cm,  L*7.14  cm 
and  the  mirror  radius  of  curvature 
Rc»500.0  cm. 


First,  we  find  solutions  to  the  wave 
equation  inside  a  perfect  infinitely-long 
rectangular  waveguide  in  terms  of  cylindrical 
waves  propagating  along  the  guide  axis.  Each 
Cartesian  component  of  the  field  must  satisfy 
the  Helmholtz  equation: 


(v7+)c7)  uix,y,z)  »  0 

Ic  ■  u/C 


(1) 


where  a  time  dependence  has  been  assinn- 

ed. 


We  are  interested  in  solutions  that 
satisfy  the  boundary  conditions  at  the  ver¬ 
tical  parallel  planes  y^O  and  ywb.  The 
proposed  solutions  are  of  the  functional 
form: 


L  B  Z 

u(x,y,z>  “  *(x,z)e  sinyy 
T  “  nn/b  ,  3  •  /)c^- 


(2) 


where  n  is  an  integer  1,2, . 

Inserting  i2)  into  ill  yields  the  following  equation  for  the  slowly  varying  complex 
function  »(x,z) 


r.he  s.cwiv  varyi.-.c  assumption  means  that  4.2^  <<  1,  and  t.herefore  we  can  neglect  - — y  in  3. 


W(2)  ■  W  /  1  +  (2/2-)  ' 

o  /  o 


BEAM  WAIST 


R(2)  ■  2(1  +  (2/2^)^) 


RADIOS  OF  CORVATURE 


2  «  3w2  RAYLEIGH  LENGTH 

o  o 

If  the  origin  of  coordinates  is  at  the  miniinuBi  spot  sice. 

In  practical  FEL  waveguide  resonator  designs,  the  transverse  dimensions  of  the  waveguide 
are  usually  much  larger  than  the  guide  wavelength.  Conseguently ,  the  amplitude  of  the  longi¬ 
tudinal  component  of  the  EjCx.y.c)  (Hjtx.y.a))  will  be  much  smaller  than  the  corresponding 
transverse  components.)^  The  assumption  of  slow  amplitude  and  phase  variations  of  v(x,2) 
with  2  in  conjunction  with  the  above  restriction  impose  on  our  theory  a  liautation  to  low- 
order  and  low- losses  inodes. 


x.i!  sin(v_y)  exp;iS_z) 
mn  n"*  n 


cos(Y„y)  exp(i3„2) 
mn  n**  n 


sin(Y  y)  exp(ia„2) 
n  sx  n-^  n 


'n'  o  X 


«x  =  =y 


H  -  -1/uu-  .“  cos(v_y)  expii3_2) 

2  o  jx  n-*  n 


Thus,  we  obtain  the  standard  paraxial  differential  equation,  the  solution  of  which  in  the 
context  of  light  beams  are  well  known-)  and  can  be  written  as 


Where : 


vx,2)  *  exp(-x^/w(2)  ^2)  exp(iS.x^/2  R(2)  -  (m*l/2)  tan" ^  ( 2/ z  )  }  H.'.x/wtzi;  N 
L  ^  Ti  o  m  m 

■  I  2®m!  /iT/l  +  Cc/Zg)  ^  j"’’  NORMALIZATION  CONSTANT 

„  ,  .  HERMITE  POLYNOMIAL 


We  break  down  our  solutions  into  two  sets  of  linearly  independent  functions:  a)  those  for 
which  Sv,  S?  <<  Ex  will  be  denoted  as  ♦mnlx>y»*>  -modes,  i. e.,  the  main  cartesian  component 
is  parallel  to  the  wide  dimension  of  the  guide;  b)  Hy,  Hj  <<Hx  will  be  designated  as 
:^(x,y,2)  in  this  case  the  dominant  electric  field  cartesian  component  is  perpendicular  to 
the  wide  dimension  of  the  guide.  The  boundary  condition  imposed  on  the  c-eigenmodes  is,  Ex, 
£z  >  0  at  y>0,  ta  for  perfect  conducting  walls;  for  9-eigenmodes  the.  boundary  condition  is 
Hv  ■  0  at  y  «  0,b.  These  arc  the  restrictions  at  the  two  parallel  horizontal  planes;  at 
t.he  vertical  planes  x  ■>  -a/2,  *a/2  the  boundary  conditions  are  automatically  satisfied 
because  all  field  amplitudes  drop  exponentially  as  a  function  of  coordinate  x  (the  scale 
being  .,(2),  the  beam  waist). 

We  present  below  a  summary  of  both  type  of  solutions. 
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Table  I. 


E  _-Modes ■« _ ex?(i8_2) 
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e£  -Modes“«,;!_exp(i3.z 
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In  "icure  I  we  illustrate  the  configuration  of  electric  field  lines  at  the  band  waist  of 
various  lew-order  modes. 


V'.N'  V 
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Figure  2.  Electric  field  linee  at  the  waist  of  low-order  waveguide  eigenaodes:  a)it-(Dodes; 
b)  s-Bodes. 


For  a  linearly  polarized  undulator  the  appropriate  modes  that  will  couple  to  the  electron 
b^am  correspond  to  mweven  and  n<^d.  However >  as  we  shall  show  in  the  next  subsection,  the 
*02]cil  ^ave  losses  that  increase  with  frequency  and  consequently  ij^hey  are  not  suitable 

for  f£l  operation  in  the  subnilliaater  region.  On  the  other  hand  the  ’‘odm*  have  a 

power  attenuation  constant  that  decreases  with  frequency. 

2.  Waveguide  losses. 

The  losses  due  to  the  finite  conductibilicy  of  the  walls  are  computed  in  the  standard 
fashion  by  evaluating  the  tangential  component  of  the  magnetic  field  H  at  the  walls;  the 
results  are, 


uc^/ugb^  (n»)* 


/ /“oS-'  - 


•  “c«/uo^ /“o*o“^  ■ 


In  Figure  3  we  plot  the  theoretical  attenuation  constant  as  given  in  Eqs.  6  as  a  function 
of  b,  the  height  of  the  resonator  channel,  and  we  cosrpare  them  to  those  of  conventional 
microwave  modes  in  rectangular  waveguides.  The  «-modes  losses  are  two  orders  of  magnitude 
below  those  of  the  more  conventional  TEQ^-mode  for  ba20  an. 

A  description  of  the  actual  mechanical 

^  ■  l.■onoo.lJl7|  design  and  paramters  of  the  UCSB  PEL 

^  - i  resonator  are  discussed  in  the  paper  at  this 

,  . -Tv  ^ — -■ —  conference  by  one  of  us  (L.E.)  and  J. 
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figure  3.  Attenuation  constant  vs.  guide 
height  b. 


Three-dimensional  numerical 


Slowly  varying  Lienard-Wiechert  solution 


The  explicit  solution  of  the  Maxwell's 
equations  for  a  point  charge  in  free  space 
can  be  written  as. 


E(t,t) 


nA{ (n-3) A6> 
R(l-3.n) 


these  are  the  Lienard-Wiechert  fields ,  where 
e  is  the  electron  charge,  1,1,  R  •  !r  -  r' 
are  respectiyely  t.he  instantaneous  normaliz¬ 
ed  yelocity  and  acceleration  of  the  particle 


the  distance  from  the  positron  of  the  particle  at  time  t'  to  the  observation  point:  and 
H  ■  R/R.  The  square  bracicet  has  the  meaning  that  the  quantity  inside  it  has  to  be  evaluated 
at  the  retarded  time  t'«t-R(t')/c. 

.  ;i  _ 

The  total  field  S(r,t)  •  2.  2i(r,t)  is  obtained  bv  summing  all  electrons  in  the  electron 

i-1 

beam.  It  will  contain  a  fast  varying  factor  of  the  form  where  u  ^  (center 

frequency  of  the  emitted  radiation)  and  uq  ■  IRc/Iq  is  the  frequency  corresponding  to  the 
period  of  the  wiggler. 


To  extract  the  slowly  varying  (slow)  part  from  E(r,t)  we  proceed  in  the  usual  fashion  of 
multiplying  by  the  fast  factor  and  then  averaging  over  many  wavelength  of  light.  This 
procedure  assumes  the  laser  is  well  above  threshold  and  that  nearly  coherence  and  monochroma¬ 
ticity  have  been  achieved. 


The  fast  part  of  the  radiation  field  can  be  visualized  as  produced  by  electrons  in  the 
combined  field  of  an  undulator  and  the  optical  field  Ej(r,t)  following  the  zero  order  tra¬ 
jectory. 


‘0=  3 


*  c3_f  e,  -  I- 


sinc)Cgt'  ei 


(8) 


where  K  »  eB^/mc^kQ,  c3q  is  the  initial  longitudinal  velocity,  is  the  undulator  wavevector 
and  3  »  Sq  sinugtej  is  the  linearly  polarized  undulator  magnetic  field.  Any  departure  from 
zero-order  trajectory,  due  to  the  ponderomotive  potential  which  force  the  electrons  to  bunch 
and  consequently  to  change  their  Sj  velocity,  will  provide  the  slow  part  cf  Che  field. 

2»T 

The  resulting  field  after  averaging  over  the  total  time  of  interaction  NT  «  —  N  (N  is  the 
number  of  periods  in  the  undulator)  is,  “ 


f(r,t)  ■  <^2(r,t)^ 


where  the  integral  is  over  the  retarded  time  t'  (time  of  evolution  for  t.he  electrons) ,  the 
sum  IS  overall  the  electrons  that  at  c«0  where  inside  of  a  wavelength  of  light  v  of  the  beam 
(Ng) ;  and  To  IS  the  retarded  time  Tq  “  /*  dt/(l-5'5). 

3 

To  derive  Eq.  9  we  assumed  electrons  pulses  nP“^**  »  NT  and  small  signal  gain  as  we  main¬ 
tained  t.he  field  experienced  by  the  electrons  constant  over  one  pass. 

To  compute  t.his  integral  we  add  the  single-particle  dynamics  described  by  the  pendulum 
equation, 
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Expanding  R  -  R^  -  ".^,5  2  and  introducing  it  in  Eq.  9  yields. 
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terms  procorticnal  to  1. RJ  are  neglected;  R^  is  the  distance  from  the  observation 
to  t.he  zero  order  oosition  of  the  oarticie  and  we  define 


wnere 


Equations  10  and  11  completely  describe  our  model^^  and  the  simultaneous  ntimerical  solution 
of  them  give  the  radiation  field  at  the  observation  point  ?  at  time  t.  The  geometrical 
arrangement  is  in  Figure  4 . 


Figure  4.  Schematic  drawing  of  a  waveguide  amplifier.  The  field  at  the  observation  point 

P  is  the  resultant  of  superposition  of  direct  wave  ard  reflections  from  the  walls. 


2.  Image  current  approach. 

Our  previous  discussion  assumed  no  spatial  constraints,  Exlr,t)  Li.  the  field  at  any  point 
r.  In^the  presence  of  perfect  conducting  walls  the  fields  must  satisfy  boundary  conditions, 
i.e.,  Etang  *  <3-  i^en  the  guide  is  rectangular  an  econctical  method  to  treat  this  ’real’ 
problem  with  boundaries  is  to  replace  it  by  another  in  an  enlarger  region  with  ’Image 
electron  beams’  and  no  boundaries. 

• 

Proceeding  as  in  electrostatic,  for  a  waveguide  of  rectangular  cross  section  v>e  ol.tain 
infinite  array  of  alternating  sign  currents.  The  total  field  at  the  observation  point  P  is, 


where  Niioage  is  the  number  of  image  currents  that  can  interfere  at  that  point.  The  radiation 
field  of  each  image  current  is  confined  to  a  cone  of  angle 
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0.75 


where  3p  •  cos~‘ (s/3^«s) ;  at  a  distance  z  the  transverse  dimensions  of  the  radiated  area  is 
ymui  "  Hence,  the  maximum  number  of  images  will  be 


•  ymax/b  ■ 


(12) 


The  UCSB  FEL  resonator  has  dimensions  a~13.34  cm  and  b>1.9  cm.  Using  the  condition  in 
Eq.  12  we  can  see  that  the  contribution  of  image  currents  with  response  to  the  vertical 
planes  are  negligible  and  it  is  sufficient  to  include  no  more  than  5  currents  in  the  (y,z) 
plane,  as  shown  in  Figure  5. 
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Figure  S.  The  problem  of  an  electron  beam 

inside  a  waveguide  is  replaced  by 
a  free-space  problem  with  'image 
electron  beams*  which  simulate  the 
conducting  surfaces. 


In  Figure  6  we  show  the  angular  distribu¬ 
tion  of  the  radiation  field  along  both  x-  and 
y-axis.  For  comparison  we  also  plot  the  free- 
space  distribution.  This  graph  shows  the 
increase  of  the  field  on  axis  respect  to  the 
free-space  case;  the  field  along  the  y-axis 
vanishes  beyond  x*1.0  cm  and  the  guided 
field  is  more  concentrated  around  the  z-axis. 

3.  Expansion  of  radiation  field  in  terms  of 
eioenmodes 


If  the  resonator  is  designed  to^operate  in 
one  of  the  transverse  eigenmode  *nn(x.y,z), 
the  presence  of  the  electron  beam,  i.e.,  the 
stimulated  radiation  produced  by  it,  will 
tend  to  create  a  steady-state  with  a  trans¬ 
verse  structure  different  from  the  initial 
one.  This  state  can  be  expressed  as  a  super¬ 
position  of  transverse  eigenmodes  of  the 
*cold*  resonator, 


mn 


Hotice  that  according  to  the  geometry  shown 
in  Figure  4  the  electrons  will  only  couple  to 
the  n-modes;  in  what  follows  we  will  suppress 
the  subscript  n. 

The  expansion  coefficients  can  be  cal¬ 
culated  in  a  semi-analytical  fashion. In 
Eq.  11  we  replace  the  free-space  Green 
function 

.  -  -iJt* 

Gj(r,r')  -  -5— 

by  the  one  appropriate  to  our  boundary  value 
problem: 
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Figure  €.  Comparison  of  the  guided  radiation 
field  in  the  x-  and  y-  directions 

with  the  free-space  case.  I,»5.7  m  y  — - - 

and  2-0.78  m.  o  -  /(x-x')^  +  (z-z')^ 


,11) 


BanXel  function  of  first  kind  of 
order  zero. 


Subsequently,  we  use  the  asymptotic  expansion  of  and  after  some  manipulations  we  obtain, 
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where  the  Fresnel  approximation  has  been  used. 

SuDStitution  of  Eq.  14  into  11  and  using  completeness  yields, 
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where 


•*  -  -i(l+2m)  +  n,c  +t—  (1-n,)  ( 

A  Zn  *  C  Z  ^ 


•  Hermite  polynomial,  nj  »  amd  Q  -  IX/(cMg);  we  have  also  assumed  a  filamentary  elec¬ 
tron  beam  on  axis. 

It  cam  be  shown  that  Agyj  decreases  as  the  Gaussian  index  m  increases. 

4 ■  Gain  per  Mode 


We  define  gain  per  mode  of  our  single-pass  amplifier  problem  as, 

i—  frif.  +2  -  lE.  I  da  j  5*  da 

G  «  ^  il-  rad'  in  ^  J  ,  I  rad  in 
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where  the  integral  is  over  the  waveguide  cross  section. 

Assuming  the  incident  optical  field  ^In  to  be  the  eigenmode  with  amplitude  on  axrs  Eg 
at  the  beam  waist,  then 

Asm 

Sn  "  2R«  - ■. - ,  :16) 

=0  l^i^mnl^ 

in*even:  n»odd. 

In  Figure  7  we  plot  the  gain  curve  of  the  three- lowest  order  Gaussian  modes  (s^0,2,41  as 
a  function  of  the  resonance  parameter  the  odd- index  modes  do  not  couple  to  the  electron 

beam  because  they  vanish  on  eocis.  Similar  arguments  show  that  the  n  index  must  be  odd. 
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Figure  7.  Gain  curves  for  different  values 
of  order  m«0,2,4  (Gaussian  order) 
with  fixed  order  n>l.  Note  that 
maximum  value  of  the  gain  curve 
decreases  and  shifts  to  a  larger 
resonance  pareuneter. 


Figure  8 . 
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Gmn  resonance  parameter.  Gam 
curves  for  different  values  of 
the  order  n«l  and  9  (vertical 
mode)  with  fixed  order  m»0 .  Note 
the  curves  are  shifted  by  an 
amount  proportional  to  (Sn-Sn-l* 
L(3-aj,es>- 


Notice  the  following  characteristics:  a)  The  curves  are  relatively  shifted  because  of  the 
phase  shifts  associated  with  each  mode.  For  the  pareuneters  (see  discussion  about  Figure  9) 
of  the  example  shown  in  Figure  7,  the  zero  of  the  m»2  mode  is  underneath  the  maximum  of  the 
iti»0  mode.  In  a  mode  competition  between  the  two,  the  m-0  will  grow  at  expense  of  the  m»2. 
bl  The  gain  of  each  mode  decreases  as  the  m- index  increases,  as  it  was  indicated  at  the  end 
of  the  previous  section.  Next,  in  Figure  8  we  varied  the  trigonometric  function  mode  index 
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In  Figure  9  w«  plot  the  value  of  the  maxi¬ 
mum  gain  as  a  f\inction  of  the  beam  waist  of 
the  initial  optical  wave.  This  result  shows 
that  there  is  an  optimal  beam  waist  at  about 
uq*1.0  cm;  this  curve  was  obtained  fixing  its 
position  at  the  center  of  the  undulator. 
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Figure  9.  Maximum  of  gain  Gwiw  vs.  ug  beam 
waist  of  input  wave.  Radiation 
field  i|  decomposed  into  eigen- 
modes  4]nn(3C,y)  which  contain  two 
parameters  ’beam  waist*  wg  and 
position  of  origin  of  coordinates 
(i.e.,  radius  of  curvature  Rvz)). 
The  beat  matching  of  radiation 
field  with  input  one  is  obtained 
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with  a  beam  waist  ug"1.0  cm. 
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